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Abstract: Peristaltic transport of a Williamson fluid 
in an asymmetric channel through porous 
medium is studied under long wavelength and 
low Reynolds number assumptions. The 
nonlinear governing equations of the peristaltic 
flow are solved using perturbation method. The 
solution for the stream function is obtained by 
neglecting inertia and curvature effects. The 
velocity distribution, the volume flow rate and the 
pressure rise are also determined. 
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I. Introduction 

Peristalsis is a well-known mechanism for 
pumping biological and industrial fluids. Even though 
it is observed in living systems for many centuries; the 
mathematical modeling of peristaltic transport has 
begun with the important works by Fung and Yih [1] 
using laboratory frame of reference and Shapiro et 
al.[2] using wave frame of reference. Many of the 
contributors to the area of peristaltic pumping have 
either followed Shapiro or Fung. Most of the studies 
on peristaltic flow deal with Newtonian fluids. The 
complex rheology of biological fluids has motivated 
investigations involving different non-Newtonian 
fluids. Peristaltic flow of nonNewtonian fluids in a tube 
was first studied by Raju & Devanathan [3]. 

Peristalsis is a mechanism adopted by many 
physiological systems and mechanical peristaltic 
pumps. Most of the physiological systems may be 
approximated as symmetric ducts. In view of this, 
Peeyush Chandra [4], Sarojamma et al [5], 
Ramachandra Rao and Usha [6], Misra & Pandey [7], 
Vajravelu et al. [8-11], Subba Reddy et al. [12, 13] 
and Srinivas et al. [14,15] made detailed studies on 
peristaltic pumping through tubes and channels. 
Brassuer and Anupampal (vide Chengel & Cimbala, 
[16]), made experiments on the mechanical 
functioning of the stomach using MRI (Magnetic 
Resonance Image). They observed that the stomach 
is a mixer, a grinder, a storage chamber, and a 
sophisticated peristaltic pump that controls the 



release of liquid and solid gastric content into the 
small intestines where nutrient uptake occurs. The 
MRI image of the stomach has revealed its 
asymmetry nature. Another physiological system 
namely uterus is also modeled as an asymmetric 
channel by Eytan and Elad [17]. These facts will 
explain the necessity of considering the physiological 
system to be asymmetric ducts also. Motivated by 
these facts, it will be interesting to study the peristaltic 
transport of Williamson fluid through an asymmetric 
channel filled with porous material. 

In this paper peristaltic pumping of Williamson 
fluid through a porous medium in an asymmetric 
channel with flexible walls is investigated. Using 
the wave frame of analysis, boundary value 
problem is solved and the results are discussed 
through graphs. 



II. Mathematical Formulation 

Let us consider the peristaltic transport of an 
incompressible Williamson fluid in a two dimensional 

channel of width d x + d 2 . The flow is generated by 

sinusoidal wave trains propagating with constant 
speed c along the channel walls. The geometry of the 
wall surfaces are defined as 



Y=H l = d l +3, cos 
Y=H 2 =—d 2 -b[ cos 

(1) 



A 



(upper wall), 



(lower wall), 



where ajand t\ are the amplitudes of the 

waves, X is the wave length, d x + d 2 is the width of 

the channel, c is the velocity of propagation, t is the 

time and .Yis the direction of wave propagation. The 
phase difference <f> varies in the range 0 < (Z) < tt in 

which <j> — 0 corresponds to symmetric channel with 

waves out of phase and ^) — n the waves are in 
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phase, further a l ,b l ,d l , d 2 and (f) satisfies the 
condition a 2 \ + b 2 i + 2~a x t\ cos<p< (d x + d 2 ) 2 . 
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Figure 1. Physical Model 

Introducing a wave frame j) moving with velocity 

c away from the fixed frame (X, Y) by the 
transformation, 

x=X-ct ,y=Y ,u=U-c ,v=V and P($=P(Xt). 

(2) 

and defining 

x y u v c~ ~h ~h X_ qf _ 

,tf3 ~ ^H* 4==^ ^ r ^ % 

^ qf c c A qf 4 /4p w up w 
A qf c 

(3) 

and using the above non-dimensional quantities, the 
resulting governing equations become (Nadeem[18]), 



du dv „ 
— + — = 0 

3* 3/ 



£Re 



3a 3a 
a — \-v- 



S 3 Re 



where 



a 2 = 



dx dy 



dv dv 
u h v — 

dx dy 



<o dz dl, 



dp 

dx dx dy 



(4) 



-o>{u+l) (5) 



dp 3 2 dT *y 3 dT yy 



dy 



dx 



dy 



(6) 



and low Reynolds number, neglecting the terms of 
order 5 and higher, equations (5) and (6) take the 
form 



dp d 



v 



du 

\dyjj 



du 



dy 



-(J 2 {u + \) 



(7) 



(8) 



3* 3/ 
^ = 0 

dy 

The corresponding boundary conditions in wave 
frame of reference are given by 

u= -1 on y=h 1 (x) 

u= -1 on y=h 2 (x) (9) 

Elimination of pressure from equations (7) & (8) yields 



dp 3 



dx dy 



f 



\+We 



v 



' du y 



ydyj 



du 



J 



dy 



-G 2 {u + \) 



(10) 



The volume flow rate q in a wave frame of reference 
is given by 

AW 
q = J udy 

h 2 (x) 

(11) 

The instantaneous flow Q(x, t) in a fixed frame is 

AW AW AW 



GM = J {u+l)dy= judy+ J ldy = 

AW AW AW q+ ( hl 



■h 2 ) 



(12) 

The time average flux Q over one period T (= XI c) 
of the peristaltic wave is 

Q=j\Qdt = ^){q+h i -h 2 )dt= q + \+d (13) 



T 



III. Perturbation solution 

Since, equation (10) is non-linear; its exact solution 
may not be possible. Therefore, we expand u, P and 
q as 

u = u 0 + Weu x +0{We 2 ), P= P 0 + We^ + OiWe 2 ), 
q=q 0 + We qi +0(We 2 ), 



(14) 



dp 



Da 



Da = -?—, e= porosity and k = permeability where P= ^ x > Substituting above expressions in 



Txx = -2[l+Wey]^ 

dx 

r xy = -[l+Wey] 



7 = 



18 



2 (du~\ 2 (du 



— + 8 — 

dx dx J 



,t w = -2[\+Wey\ 



dx 

equation (10) and boundary conditions (9), we get the 
following system. 

System of order We 



^dxj 



+ 



V 



dy 



+ 8' 



dv 
dx 



+ 28 2 



J 



r d_y} 

\dyj 



(15) 



Here 5, Re, We represent the wave, 
Reynolds and Weisseing numbers, respectively. 
Under the assumptions of long wavelength 5 « 1 



dp 0 = d 2 u { 
dx dy 2 

and the respective boundary conditions are 

u 0 =-l for y = h l5 u 0 =-1 fory = h 2 

(16) 
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System of order We 1 

dp x _ d 2 u x d (da ^ 



+ -a L u x (17) 

dx dy dy{ dy ) 

Uj = 0, fory^hj ,1^=0, fory = h 2 (18) 

T T T 0 

Solution for system of order We 

Solution of Eq. (15) satisfying the boundary 
conditions (16) can be written as 



(19) 
where 



a 



f f 
v 



py±ir<r°y 



e" J +ke " J -\ 



J 



-1 



crtu ah 
e -e 1 



J 



a{h 2 -h^) (j{\-h 2 ) 



- e 



, and the volume flow 



J 



rate q 0 is given by 



f 



-T- 



From Equation (20), we get 



(20) 



\-ke° h ) ( e° k -e ah ) -Jt( e ak -e° h ) -<r(4 -A, ) 



(21) 

Solution for system of order Wfe 1 

Substituting the zeroth-order solution (19) into 
(17), the solution of the resulting problem satisfying 
the boundary conditions take the following form. 



( A+B+ Qe^+I^E+F-Gje 



-Fe^e 2c7J/ + 



3a 3 



D 

v 



D 



1 



o 2 



J 



(22) 



where 



A =E(e^-e^),B 



C= 2 (dp/ 



K dx j 



3o 



k 2 ^<7(A 2 -2;0 _ e< r(A l -2A 2 ) j 



D = e^" A) -e (j(A ^ ) ,E=^r^ 

V dx) 



a 



F= 2 



3a 3 



(l-i^)'. G= ^3 



2 r^v 



V (fry 



2 -2(7/^ 



re 



and the volume flow rate q± is given by 



A 

a = — e 
Hl D 



-2ot\ 



+Ee 



' e — e 1 
v ° J 



A 
+ — e 
D 



-2ol\ 



V 



a 



-E{l\-h 2 ) 



J 



+H 



(23) 



where 

H= 



(B+Q 



D 









(n-q 




I <* ) 


I ° J 


D 


[ <* J 



F 



dx3cr 



From equation (23), we get 



where 



d Pl _a 3 ( qi -H)D 



dx 



I 



(24) 



Substituting equations (21) and (24) in to equation 
(14) and using the relation (14), we get 



ck~ J 
25) 

T I 
where J = — , 

D 

one wavelength, we get 

(21) 

idc „ J J 



J) 



( 



Integrate above equation over 



qih-k))- <fm 4+if+4-m+Hh-}iH 



(26) 

where 1=^ + 4 + 4, - 4, + 4 + 4- ^- 4 



4 = ^-(i- J te-"M ! 



2 ^ -2ahj a(2h2-hj) ah2 ^ 



-e 
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4=^(i-^)V 2 *v 
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Fig. 2. The variation of AP with Q for different values of (T 
for a fixed a = 0.8, b =0.5, (f) =0, d =2; We = 0.03 



IV. Results and Discussions 

From equation (26) we have calculated the 
pressure difference as a function of Q for different 
values of permeability parameter a and different 
phase differences (f) for a fixed a=0.8, b=0.5, d=2, 

We=0.03 and is shown in figures (2) and (3).We 
observe that the larger the parameter a the greater 
the pressure rise against which the pump works. We 

observe that for a given AP , the flux Q increases 
with increasing a .For free pumping there is no 

difference in flux Q for increase in (j.We also 
observed that the pressure rise increases with 
increasing phase difference <p . 

The variation of pressure rise with time 
averaged flow rate is calculated from equation (26) for 
different values of We and different phase differences 
for a fixed a=0.8,b=0.5,d=2, <T=1.5 and is shown in 
figures (4) to (6). We observe that the larger the We, 
the smaller the pressure rise against which the pump 
works. We observe that for a given AP,the flux Q 

decreases with increasingWe. For a given flux^ the 

pressure rise decreases with increase We. 

The variation of AF with time averaged flow 
rate is calculated from equation (26) for different 

values of the phase difference 0, for a fixed a=0.8, 

b=0.5,d=2, (j=0.5, We=0.03 and is shown in figure 
(7). We observe that the larger the phase 
difference the greater the pressure rise against 

which the pump works. We observe that for a 
given AP, the flux Q increases with increasing 0. 

For a given flux Q the pressure rise increases with 
increasing 0. 
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Fig. 3. The variation of AP with q for different values of <T 

for a fixed a = 0.8,b =0.5, <p = 71/6, d =2; We = 0.0 




Fig. 4. The variation of AP with Q for different values of We 
for a fixed a = 0.8,b =0.5, (f> = 0, d =2; G = 1 .5. 





Fig. 5. The variation of AP with Q for different values of We 
for a fixed a = 0.8, b =0.5, </> = Til 6 , d =2, G = 1 .5. 
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Fig. 6. The variation of AP with Q for different values of We 
for a fixed a = 0.8, b =0.5, (j) = n/3, d =2, G = 1.5. 
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Fig. 7. The variation of AP with Q for different values of 0 
for a fixed a = 0.8, b =0.5, We =0.03, d =2, G = 0.5. 
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